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We study the nuclear ground-state properties by using the unitary-model-operator 
approach (UMOA). Recently, the particle-basis formalism has been introduced in the 
UMOA and enables us to employ the charge-dependent nucleon-nucleon interaction. 
We evaluate the ground-state energies and charge radii of 4 He, 16 0, 40 Ca, and 56 Ni 
with the charge-dependent Bonn potential. The ground-state energy is dominated by 
the contributions from the one- and two-body cluster terms, while, for the radius, the 
one-particle-one-hole excitations are more important than the two-particle-two-hole exci¬ 
tations. The calculated results reproduce the trend of experimental data of the saturation 
property for finite nuclei. 

Subject Index D10 


Introduction. In nuclear theory, one of the most fundamental problems is to describe and 
understand the nuclear structure based on nuclear forces. This attempt has been made 
possible owing to the progress of the computational power and that of nuclear many-body 
approaches. For light nuclei with the mass number around 4 < A < 12, the ah initio methods 
such as the no-core shell model (NCSM) [1] and the Green’s function Monte Carlo method [2] 
have been applied extensively to the ground- and excited-state properties. However, it is 
difficult to apply these methods to heavier nuclei because of the current computational lim¬ 
itations. In order to investigate the medium-mass region, one can apply the many-body 
methods such as the the coupled-cluster method (CCM) [3] and the unitary-model-operator 
approach (UMOA) [4]. By employing these methods, ground-state energies of doubly magic 
nuclei have been calculated [5-10]. In the CCM, the systematic calculations have been per¬ 
formed for the ground- and excited-state energies around oxygen [11] and calcium [12]. Such 
studies have been also done in the self-consistent Green’s function method [13,14] and the 
in-medium similarity renormalization group approach [15-17]. 

Besides the nuclear many-body approaches, one of the other important ingredients of 
numerical calculations is the nuclear interaction or potential. Recently, modern nucleon- 
nucleon (NN) potentials [18-22] have been developed and can reproduce with high precision 
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the proton-proton and proton-neutron scattering phase shifts, as well as the deuteron proper¬ 
ties. It is known that these NN potentials generally have non-perturbative properties caused 
by the strong repulsive core and the singularity coming from the tensor force at short dis¬ 
tances. Therefore, it is difficult to apply such bare potentials directly to the nuclear structure 
calculations because they require the huge model space spanned by single-particle basis 
states. In order to tailor these short-range correlations to be computationally tractable, the 
effective-interaction theory is generally introduced. In the present work, we use the UMOA, 
which is one of the energy-independent Hermitian effective-interaction theories. The calcula¬ 
tions of ground-state energies and charge radii under the isospin symmetry were done for 16 O 
and 40 Caby the UMOA [4,23,24]. Recently, the particle-basis (proton-neutron-basis) formal¬ 
ism was introduced in the UMOA [5]. Since the protons and neutrons are treated explicitly 
as different particles in the particle-basis formalism, charge-dependent NN potentials can be 
used. The ground-state energies of 16 0, 40 Ca, and 56 Ni are discussed by using such a formal¬ 
ism in Ref. [6]. However, the charge radius, which is indispensable for understanding of the 
nuclear saturation property, has not been calculated with the charge-dependent framework 
in the UMOA. Nowadays, charge radii of many nuclei across the nuclear chart have been 
measured through the elastic electron scattering or the isotope shift [25]. To investigate the 
charge radius, it is, in principle, desirable to take into account the charge-dependence of the 
NN interaction by the particle-basis formalism. 

Moreover, the studies of the infinite nuclear matter system at high densities, especially 
the equation of state, are necessary to investigate the astrophysical objects such as neutron 
stars and supernova explosions. To obtain the equation of state, one of the important issues 
is to satisfy the saturation property of the infinite nuclear matter. However, it is difficult 
to reproduce microscopically the saturation property of such system, at least without use 
of genuine three-nucleon forces. In order to obtain more reliable predictions for the infinite 
nuclear matter system, one can start from the reproduction of the properties of finite nuclei, 
because the experimental data are available only for finite nuclei. The saturation property 
of 16 0 and 40 Ca was investigated with Bonn potentials in the Briickner-Hartree-Fock (BHF) 
method, and the influence of effective masses of mesons were discussed [26]. Furthermore, the 
relativistic effect on 16 O and 40 Ca was included through the Dirac-Briickner-Hartree-Fock 
method [27] and improves the BHF results. On the other hand, it is shown that the results for 
such nuclei in the UMOA are closer to the experimental values than the BHF results [4,24] 
even without relativistic effect. Before discussing quantitatively the various contributions, 
for instance, from the relativistic effects and three-body forces, in this study, we employ 
the CD-Bonn potential [20], which is one of the modern high-precision NN interactions and 
holds the charge-symmetry and charge-independence breakings. In this Letter, we calculate 
the ground-state energies and charge radii, corresponding to the inverse of density, of 4 He, 
16 0, 40 Ca, and 56 Ni, and discuss the saturation property of finite nuclei. 


Framework. Here, we briefly present the calculation method of the ground-state energy 
and charge radius in the UMOA with the particle-basis formalism. The details can be found, 
for example, in Refs. [4, 23, 28]. The essence of the UMOA is to construct the effective 
Hamiltonians from the original total Hamiltonian H and to obtain the physical observables 
with effective operators. The transformed total Hamiltonian, H = e~ s He s , is constructed 
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by the unitary operator e s , and takes into account two-particle-two-hole (2p2h) excita¬ 
tions considered as the most important correlations in the nuclear structure. The exponent 
S is called as the two-body correlation operator. Since e s is the many-body operator, H 
becomes the many-body operator. The transformed Hamiltonian can be decomposed into 
H = iT 1 ) + H^ + H ® + • • •, according to the number of interacting particles. Here, H^ l \ 
H (2 ), and are the one-, two-, and three-body cluster terms, respectively. In the prac¬ 
tical calculation, we treat the one- and two-body cluster terms directly and obtain S by 
solving the decoupling equation, QH^P = 0. Here, P and Q are the projection operators 
onto the OpOh and 2p2h states, respectively. Note that the effective interaction is determined 
self-consistently with the one-body potential. 

To take the other important correlations into account, we diagonalize the transformed 
Hamiltonian in the OpOh and lplh space. This procedure is needed because the one-body 
cluster term induces the lplh excitations. After the diagonalization, we have the energy 
of the one- and two-body cluster terms, E 1+2BC . Due to the diagonalization, the ground 
state is expressed as the linear combination of the unitary transformed wave functions in the 
OpOh and lplh space, [T = e s (go + J2 a /3 dapOabp )|</>o • Here, |</>o is the particle-hole vacuum 
state, and go and g a p are coefficients of the OpOh (particle-hole vacuum) and lplh states, 
respectively. Subscripts a, /3, • • • denote the set of harmonic-oscillator quantum numbers. The 
Oq and b ^ are creation operators of the particle state a and the hole state /3, respectively. 
In addition to E l+2BC , we evaluate the contribution of the three-body cluster term to 
the ground-state energy, E 3BC , through order S 2 [4,23,28] and obtain the total energy, 
E„, s , = E 1+2BC + E 3BC . 

The charge radius R c is evaluated by using the equation, R 2 = (’F|r 2 |'F) + R 2 + ( N/Z)R 2 
[29]. Here, R 2 and R 2 are the squared radius obtained from the experimental charge dis¬ 
tribution of the proton and neutron, respectively. Here, we take i? 2 = 0.832 fm 2 [30] and 
R 2 = —0.115 fm 2 [25]. The Z and N are the proton and neutron number, respectively. 
What we calculate in the UMOA is ('k|r 2 |T) defined by ~ Rc.m) 2 \^)/Z. The 

Ti and R c .m. are the coordinate vectors of the ith proton and the center-of-mass (c.m.) 
of A nucleons, respectively. For N = Z nuclei, R c . m . is approximately equal to the c.m. 
coordinate vector of protons, whose expectation value is approximated as Os state, i.e., 
(\k|i? 2 m ['I') ~ 3h/(2Amuj). Then, the term ^T(T|r* 2 |T)/Z is expanded up to order S 2 into 

Z 2 2 

(^1 ^ ^ r j 1^) — 90 ^ ^ r Oi , a' H ^ ^ S a p a 'f3'S^p a '(3'V a ry ~ ^ ^ S a p a 'p'S a pyp'T a 'y 

i =1 a'<pF a/3 r y>pF a(3>pF 

a'f3'<pF cZ (3 ,r y'<pF 

^90 ^ ^ Qaa'^aa' ^90 ^ ^ ^ol(3ol' f3'^aa' f zpr ( 1 ) 

ol>pf ol(3>pf 

a'<pF a'f3'<pF 

H - ^ ^ QolOl' QolOl' T/3'ft' E 9aa' 9^a'^af3 ^ ^ 9olol' Qa/S'^a.'(3' • 

a>pF ap>pF ol>pf 

OL'j3'<pF ol'<Pf ol'P'<pf 


Here, we use the notations, r a p = {a\r 2 \f3) and S a /3^s = {&/3\S\'y5). The pf denotes the Fermi 
level. The matrix element of r 2 is evaluated only with respect to the proton states. Note 
that summations in Eq. (1) run over only the proton states. However, indices j3 and (3 1 in 
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Table 1 Ground-state energy for each nucleus. We take p\ = 12 and fruj m i n = 18 MeV 
for 4 He, pi = 14 and huj m i n = 15 MeV for 16 0, p\ = 18 and frw m i n = 14 MeV for 40 Ca, and 
pi = 20 and huj m i n = 14 MeV for 56 Ni. The definitions of E 1+2BC , E 3BC , and E gs . are given 
in the text. The experimental values are taken from Ref. [34], All the energies are in units 
of MeV. 



4 He 

16 0 

40 Ca 

56 Ni 

^1+2 BC 

-26.13 

-115.58 

-334.36 

-454.84 

E 3BC 

-1.60 

-3.82 

-5.92 

-18.20 

Eg.s. 

-27.73 

-119.39 

-340.28 

-473.04 

Expt. 

-28.30 

-127.62 

-342.05 

-483.99 


the second, third, and fifth terms in Eq. (1) run over both of the proton and neutron states. 
The isospin factor f Za , is —1 (+1) for proton (neutron). 

Results and discussions. In this section, we show the results of the ground-state energy 
and charge radius of 4 He, 16 0, 40 Ca, and 56 Ni, and discuss the saturation property of these 
nuclei. All the calculated results are obtained with the CD-Bonn potential [20]. The Coulomb 
interaction is included in the proton-proton channel. Following Refs. [5,6], we adopt the two- 
step decoupling method. At the second-step decoupling, we add the term, i7 cm . = /3(T c . m . + 
17c. m . — 3hw/2), to the intrinsic Hamiltonian H — T c . m . so as to remove the c.m. spurious 
motion by the Glockner-Lawson prescription with [3 = 1. Here, T c . m . and t/ c . m . are the kinetic 
and harmonic oscillator potential of c.m., respectively. The details can be found in Ref. [5]. 
The model-space size is defined by the sum of the relevant quantum numbers of the two-body 
state, pi = 2 n a + l a -f 2 np + Ip. Here, n a and l a are the principal and azimuthal quantum 
numbers of the harmonic-oscillator state a, respectively. We perform the calculations with 
various pi and the harmonic-oscillator energy hu, and investigate the pi- and hu-dependence 
of the ground-state energies and charge radii. 

We observed that the ground-state energy lowers monotonically as increasing pi similar to 
Ref. [6]. This is non-trivial because our calculations do not have to preserve the variational 
principle mainly due to the truncation of four- and higher-body cluster terms. When many- 
body cluster terms are truncated, the dependence of the model-space size on the energy does 
not obey the variational principle and can be found, for example, in the case of the NCSM [31] 
with the Lee-Suzuki transformed effective interaction [32,33]. Actually, our ground-state 
energy of 4 He is overbound a little compared with the NCSM result in the sufficiently large 
model space, which is thought of as the exact solution with the CD-Bonn potential. Our 
final results of the ground-state energy depend slightly on frw even if the pi-dependence 
vanishes. Since the choice of Huj is arbitrary in nature, the htc-dependence should vanish 
at sufficiently large p \. In the case of the CCM, the hw-dependence at the coupled-cluster 
double (CCD) level, similar to our framework, vanishes at the coupled-cluster single and 
double (CCSD) level [35]. Therefore, the hic-dependence in the UMOA could be reduced 
by introducing the one-body correlation operator, in addition to the two-body correlation 
operator, it was also shown that the ground-state energy at fku m ; n in the CCD is close to 
the result in the CCSD. Here, huj m ; n is the value of ftw minimizing the ground-state energy. 
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Table 2 Charge radius R c for each nucleus. The entries of “ R c w/o lplh and 2p2h 
correlations”, U R C w/o lplh correlations”, and ll R c w/o 2p2h correlations” are the results 
from Eq. (1) with both of S a p a 'p' = 0 and g aOL ' = 0, g aa > = 0, and S a p a 'p' = 0, respectively. 
The model-space size and hu m ; n are same as in Table 1. The experimental values are taken 
from Ref. [25]. All the radii are in units of fm. 



4 He 

16 0 

40 Ca 

56 Ni 

R c w/o lplh and 2p2h correlations 

1.81 

2.59 

3.08 

3.28 

R c w/o lplh correlations 

1.82 

2.60 

3.09 

3.29 

R c w/o 2p2h correlations 

1.67 

2.44 

2.97 

3.20 

Rc 

1.67 

2.44 

2.97 

3.19 

Expt. 

1.68 

2.69 

3.48 



Thus, we tabulate the ground-state energy of each nucleus at Hw m { n in Table 1. The E 1+2BC 
is the energy obtained from the one- and two-body cluster terms. The contribution of the 
three-body cluster term, E iBC , is much smaller than e 1+2BC and attractive for all nuclei 
examined here. This tendency is also observed in the calculated results obtained by the 
CCM [8,9]. The contributions from the higher-body cluster terms would be less than those 
of the three-body cluster term. Thus the ground-state energies, Fg.s. = e 1+2BC + E 3BC , are 
expected to almost converge with respect to the cluster expansion. 

In contrast to the ground-state energies, charge radii have the sizable ^-dependence 
except for 4 He even if the results converge with respect to p \. The charge radii decrease 
monotonically as Hu increases. For example, the charge radius of 16 0 ( 56 Ni) changes from 
2.76 (3.56) fm at Hu = 11 MeV to 2.27 (2.87) fm at Hu = 20 MeV. To evaluate the charge 
radius in the UMOA, Hu was formerly taken to Hu m i n [4,24], Recently, it was demonstrated 
that the radius at Hu m \ n coincides the ^-independent result in the CCM [35]. Therefore, 
in Table 2, the charge radius of each nucleus at Hu m i n is tabulated. The contributions of 
the correlation operator S and the diagonalization coefficients g aa > are also shown as those 
of 2p2h and lplh correlations, respectively. Since the effects of the correlation operator 
S are much smaller than those of g aa ', lplh excitations to the charge radius are more 
important than 2p2h excitations. In other words, for the radius, it can be expected that the 
consideration of the one-body correlation operator is more important than that of the two- 
body correlation operator. All the charge radii investigated here shrink compared with the 
experimental values as increasing the mass number. Note that the charge radius of 56 Ni has 
not been measured yet. In our framework, neutron and matter radii can be also calculated in 
the same way. The differences between the proton and neutron radii for the nuclei examined 
here are smaller than 0.02 fm. Therefore, we only discuss the proton radius in this work. 

In Fig. 1, the saturation property of 4 He, 16 0, 40 Ca, and 56 Ni is illustrated. Our calculated 
results of 16 0 and 40 Ca are consistent with the results by the earlier UMOA calculations, 
which were obtained with the realistic NN potentials at hu m i n [4,24], The dashed curve, 
which is obtained by the empirical formula given by Bethe and Weizsacker [36], E/A = ay + 
asA” 1 / 3 + ac^ 2 /A 4 / 3 + a\(N — Z) 2 /A 2 — apA~ 7 / 4 , and the charge radius, R c = roA 1 / 3 [36, 
37] shows the systematic behavior for light- and medium-mass N = Z nuclei. Here, we adjust 
ay = —16.10 MeV and ro = 1.05 fm, and use as = 18.56 MeV, ac = 0.717 MeV, a\ = 28.1 
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MeV, and ap = 34.0 MeV as found in Ref. [36]. Our calculated results show similar pattern 
to this empirical curve. However, taking a closer look at our results with the NN interaction 
only, all the binding energies are underbound to and also all the charge radii are smaller 
than the experimental data. When we use the other realistic NN potentials, the binding 
energies are expected to increase (decrease) with decreasing (increasing) charge radii [26], 
and distribute over the Coester line [36]. Actually, the numerical results by the UMOA 
given in Refs. [4,24], show such a correlation between the binding energy and charge radius. 
Therefore, it is likely that the results do not approach the experimental data even if we 
use other modern high-precision NN interactions. The deviation from the experimental data 
might come from the lack of the genuine three-nucleon force effect, because it gives the 
attraction in light nuclei as shown in the Green’s function Monte Carlo method [38]. Also, the 
three-nucleon force through the relativistic framework stretches the nuclei as found from the 
comparison between the Briickner-Hartree-Fock and Dirac-Briickner-Hartree-Fock results 
[27]. Combining these facts, the three-nucleon force seems to be necessary to reproduce the 
binding energy and charge radius simultaneously. From the recent ab initio calculation of 
Green’s function theory [39], the three-nucleon force effect actually increases the matter radii 
of 16 0 and 44 Ca. However, the discrepancy between the experimental data and the recent ab 
initio results still remains and needs to be further investigations [14,40]. As for 4 He, we can 
compare with the result obtained by the NCSM [31] which is also plotted in Fig 1. Our result 
is close to the NCSM result. The NCSM result can be considered as the exact solution with 
the CD-Bonn potential, because the calculation was performed at sufficiently large model 
space where the fvjj- and model-space size dependencies are negligible. Our charge radius for 
4 He can be determined uniquely, because of the weak foj-dependence. For 16 0, 40 Ca, and 
56 Ni, charge radii have sizable ^-dependence. Thus, we cannot discuss these results in the 
same accuracy as 4 He. It is necessary to obtain the practically Iku-independent results. To 
achieve this, the consideration of the one-body correlation operator would be desirable, i.e. 
the construction of the effective Hamiltonian which does not induce the lplh excitations in 
addition to the 2p2h excitations by the unitary transformation. 

Summary. We represent the way to obtain the ground-state energy and charge radius in 
the particle-basis formalism of the UMOA. We evaluate the ground-state energy and charge 
radius for 4 He, 16 0, 40 Ca, and 56 Ni with the CD-Bonn potential. The results are obtained 
in the model space defined by the two-body state, p\ = 2 n a + l a + 2 np + Ip, and converge 
with respect to p\. We find almost converged ground-state energies with respect to the 
cluster expansion, because the contribution of the three-body cluster term estimated here 
is much smaller than the one- and two-body cluster terms. All the ground-state energies 
calculated here are underbound to the experimental data. The contribution of the two-body 
correlation operator S to the charge radius is less important than lplh excitations. All the 
charge radii estimated here are smaller than the experimental data. From our results of 
the UMOA, we show the obtained saturation property of finite nuclei, consistently with the 
trend of experimental data. Our result of 4 He is close to the ab initio solution with the same 
interaction. 

In this Letter, we take frw minimizing the ground-state energy. The validity of this choice 
is discussed in the CCM [35]. However, it is necessary to obtain the virtually ^-independent 
results for the quantitative comparison with the experimental data. In the CCM, the role 
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Fig. 1 Ground-state energy per nucleon and inversed charge radius of 4 He, 16 0, 40 Ca, 
and 56 Ni. Filled symbols correspond to the experimental values. As for the charge radius of 
56 Ni , we substitute the experimental value of 58 Ni. Open symbols correspond to the numerical 
results by the UMOA. The cross symbol is the NCSM result [31]. The circles, down triangles, 
up triangles, and squares are for 4 He, 16 O, 40 Ca, and 56 Ni, respectively. The dashed curve 
is obtained from the Bethe-Weizsacker formula for the energy and R c = roA -1 / 3 for the 
radius [36,37]. The adopted parameters are shown in the text. 


of the one-body cluster operator is investigated [35], and it is shown to reduce the /in¬ 
dependence. If we include the one-body correlation operator, in addition to the two-body 
correlation operator, the Au-dependence is expected to be weakened. The results in the 
UMOA with the one-body correlation operator will be reported elsewhere in the near future. 
It is also found that the genuine three-nucleon force seems to be necessary so as to quanti¬ 
tatively reproduce the saturation property of finite nuclei. The inclusion of the three-body 
force in the UMOA is under way. 
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